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We present a search for decays of B mesons to final states with a b\ meson and a p or if* (892) 
meson. The search is based on a data sample consisting of 465 million BB pairs collected by the 
BABAR detector at the SLAC National Accelerator Laboratory. We do not observe any statistically 
significant signal. The upper limits we set on the branching fractions range from 1.4 to 8.0 x 10~ 6 
at the 90% confidence level (C.L.), including systematic uncertainties. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30. Er 



Measurements of charmless hadronic B decays are a search for new physics effects. One of the outstanding 
powerful tool to test standard model predictions and problems is represented by the so called polarization puz- 
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zle in decays of B mesons to a pair of spin-one mesons. 
Simple helicity arguments predict a longitudinal polariza- 
tion f L close to 1. Contrary to this, several vector- vector 
(VV) decay modes such as B -> 4>K* Q. /i • p+ivT* 
0, and B — * wK* 3] exhibit / L - 0.5. Possible ex- 
planations for this puzzle have been proposed within the 
standard model 0] and in new physics scenarios [B|. 

The measurement of the branching fractions and po- 
larization of charmless decays of B mesons to an axial- 
vector and vector meson (AV) may shed light on the size 
of the amplitudes contributing to charmless £>-meson de- 
cays and on their helicity structure. Theoretical predic- 
tions of decay rates have been performed with the naive 
factorization (NF) 6] and QCD factorization (QCDF) 
approaches. The NF calculations find the rates of 
B — ► AV decays to be smaller than the correspond- 
ing B decays to an axial- vector and pseudo-scalar meson 
(AP). The more complete QCDF calculations find the 
reverse, primarily due to the larger decay constants (p 
vs 7r for instance); the expected branching fractions for 
the AV modes are substantial in several cases, as large 
as 33 x 1CT 6 for the B° -> b~[ p+ final state. 

Additionally, decays of B mesons to charmless AV final 
states may be sensitive to penguin annihilation effects, 
which tend to enhance certain modes while suppressing 
others. It is thus important to investigate the largest 
possible number of final states. 

Measurements of the branching fractions to AP modes 
b\h, where h denotes a charged or neutral pion or kaon, 
are presented in Ref. @. The results are in good agree- 
ment with the predictions of QCDF [§] . Searches for the 
AV decays to the final states afp T and afK*° are pre- 
sented in Ref. [Icj . with upper limits on the branching 
fractions of 30 x 10~ 6 and 1.6 x 10~ 6 (at the 90% C.L.), 
respectively. In this paper we search for all charge com- 
binations of decays of a B meson to a final state contain- 
ing a bi meson and a p or if* (892) meson. No previous 
searches for these decays have been reported. 

The data sample used for these measurements was 
collected with the BABAR detector at the PEP-II asym- 
metric e + e~ collider located at the SLAC National Ac- 
celerator Laboratory. The integrated luminosity taken 
at the T(AS) resonance (center-of-mass energy y/s = 
10.58 GeV) corresponds to 424 fb _1 and is equivalent 
to (465 ± 5) x 10 6 BB pairs. The BABAR detector is 
described in detail elsewhere (l2] |. 

We reconstruct B-meson daughter candidates through 
the decays b\ — ► coir (we assume this branching frac- 
tion tO be 1 [IH), CJ -> TT+TT-TT , p+ -> 7T+7T , p° -> 

7T+7T-, K*° -► K+tt-, and K*+ -> K+n° or K°tt+. 
We impose the following requirements on the masses of 
the selected candidates: 1000 < m(bi) < 1550 MeV, 
740 < m(uj) < 820 MeV, 470 < m{p) < 1070 MeV, 
and 755 < m(K*) < 1035 MeV; these cuts allow some 
sidebands, which help estimating the background level. 
Neutral pions are reconstructed via the decay 7r° — > 77; 



photon candidates with a minimum energy of 50 MeV are 
combined, and we require the pion energy to exceed 250 
MeV in the laboratory frame. The invariant mass of the 
7T° candidate is required to be in the interval 120—150 
MeV. We select Kg — > ir + ir~ candidates in the mass 
range 486 < m(K®) < 510 MeV; a kinematic fit con- 
straining the two pion tracks to originate from the same 
vertex is performed and we require the K® flight length to 
be greater than three times its uncertainty. The daugh- 
ters of bi, u>, p and K* are rejected if their particle iden- 
tification signatures are consistent with those of protons 
or electrons. K + candidates must be positively identified 
as kaons, while 7r + must fail kaon identification. Unless 
otherwise stated, charge-conjugate reactions are implied. 

The helicity angles of the (axial-) vector mesons are 
measured in their rest frame. For the b\ candidate, the 
helicity angle is defined as the angle between the flight 
direction of the pion from the b\ — ► lutt decay and the 
direction of the boost to the b\ rest frame. We define 
the helicity angles of the p and K* mesons in an anal- 
ogous manner using the direction of the daughter pions 
[for the p ± (p°) we use the (positively) charged pion]. 
Finally, the helicity angle of the w is taken as the an- 
gle between the normal to the 3tt decay plane and the 
direction of the boost to the u) rest frame. To suppress 
backgrounds originating from low- momentum particles, 
we apply the selection criteria summarized in Table [J 
Integration over the angle between the b\ and V decay 
planes yields the following expression for the distribution 
F(9a,0 v ) cx cPT I d cos 6 Ad cos 6y in the b\ and p/K* he- 
licity angles 9a and 9y\ 



F(8 Al 9 v ) = f L 
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2 


cos 2 9 a + 


C Q 


sin 2 9a 



(Wl) 







sin 2 9 a + 


C 





(1 + cos 2 9 A ) 



COS 9y + 



sin 2 (9y. (1) 



Here is the longitudinal polarization fraction 

I j4q 1 2 / \Ai\ 2 , where Ai, i = —1, 0, 1, is a helicity ampli- 
tude of the B — > AV decay. The d are the helicity am- 
plitudes of b\ — » urn; by parity conservation C_i = C\. 
The 61 decays have been studied in terms of the two 
parity-allowed S and D partial wave amplitudes, which 
have the measured ratio D/S = 0.277 ±0.027 [11]. From 
this we obtain the ratio of helicity amplitudes in Eq. Q] 

m 

Ci = l + (£>/S)A/2 
Co 1 - V2(D/S) ' 

Two kinematic variables characterize the decay of a B 
meson: the energy-substituted mass tubs = i/s/4 — 
and the energy difference AE = Eb — where 
(Eb,Pb) is the B- meson four-momentum vector ex- 
pressed in the T(4S) rest frame. The correlation between 
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TABLE I: Selection requirements on the helicity angles of 
B-daughter resonances. 



State 


p/K* helicity 


6i 


helicity 


b-p+ 


-0.50 < cos 


6 P < 1.00 


-1.0 < 


COS 


< 1.0 


6?p+ 


-0.50 < cos 


P < 0.80 


-1.0 < 


cos 8t 1 


< 0.6 


bt P ° 


0.0 < |cos( 


1 p \ < 0.85 


-1.0 < 


cos 9b 1 


< 1.0 


6?P° 


0.0 < cos( 


? p < 0.85 


-1.0 < 


cos 6b 1 


< 0.7 


bfK* 
blK* 


-0.85 < cos( 
-0.85 < cos ( 


?jc- < 1.0 

?K» < 1.0 


-1.0 < 
-1.0 < 


cos 

COS Sfcj 


< 1.0 

< 0.8 



the two variables is at the few percent level. The reso- 
lution on ibes is about 2.6 MeV, while the resolution 
on AE varies between 20 and 40 MeV depending on the 
number of ir° mesons in the final state. We select events 
with 5.25 < m ES < 5.29 GeV and \AE\ < 0.1 GeV ex- 
cept that for b\p + we require -0.12 < AE < 0.10 GeV 
to allow for the broader signal distribution when two ir° 
mesons are present. The average number of B candidates 
per event in the data is between 1.3 and 1.6. We choose 
the candidate with the highest value of probability in the 
fit to the B vertex. 

The dominant background originates from continuum 
e + e~ — > qq events (q = u,d,s,c). The angle 9t be- 
tween the thrust axis |14] of the B candidate in the 
T(AS) rest frame and that of the remaining particles 
in the event is a powerful discriminating variable to sup- 
press this background. Continuum events peak near 1.0 
in the | cos 0t | distribution, while B decays are almost 
flat. We require | cos 9t\ < 0.7 for all the decay modes ex- 
cept p° for which we require | cos9t\ < 0.55, because 
of substantially higher backgrounds. To further reduce 
continuum background we define a Fisher discriminant 
(J 7 ) based on five variables related to the event topol- 
ogy: the polar angles, with respect to the beam axis, of 
the B candidate momentum and the B thrust axis; the 
zeroth and second angular moments Lq and L2 of the 
energy flow, excluding the B candidate; and the flavor 
tagging category The first four variables are calcu- 
lated in the T(4S I ) rest frame. The moments are defined 
by Lj = ^2iPi x |cos#j| J , where 0j is the angle with re- 
spect to the B thrust axis of track or neutral cluster i, pi 
is its momentum. The Fisher variable provides about one 
standard deviation of separation between S-decay events 
and combinatorial background. 

The signal yields are obtained from extended maxi- 
mum likelihood fits to the distribution of the data in nine 
observables: AE, tties, rrik, and cos 0k] rrik and 9k are 
the mass and the helicity angle of meson k (k = b\, u>, 
and either p or K*). For each category j (signal, qq 
background and backgrounds originating from BB de- 
cays), we define the probability density functions (PDFs) 
Vj(x) for the variable x, with the associated likelihood 



C: 




i=i j 

vAurypMvjyPiV*) II (^(K)^(cos^)) , 

k 

where Yj is the event yield for component j and N is the 
number of events entering the fit. We separately model 
correctly reconstructed signal events and self-crossfeed 
(SXF) events, which are signal events for which parti- 
cles are incorrectly assigned to the intermediate reso- 
nances, or particles from the rest of the event are se- 
lected. The fraction of SXF is 0.33-0.57 depending on 
the final state. The signal yields for the branching frac- 
tion measurements are extracted with the use of correctly 
reconstructed signal events only. 

Backgrounds originating from B decays are modeled 
from Monte Carlo (MC) simulation [16(. We select the 
most significant charmless modes (20-40 for each signal 
final state) entering our selection and build a sample tak- 
ing into account measured branching fractions or theoret- 
ical predictions. The expected charmless BB background 
yield varies between 26 and 330 events, depending on the 
final state. The samples include the nonresonant contri- 
butions affecting b\p (61 if*), measured in our data by 
fitting the central regions of the 6i7T7r (biKir) and u)irp 
(umK*) Dalitz plots. We assume the probability of the 
four-body nonresonant contributions to pass our selec- 
tions to be negligible. We do not introduce a component 
modeling B decays to charmed mesons, since this back- 
ground is effectively absorbed in the qq component. 

For the K* modes we consider the potential back- 
ground contribution originating from Kn S'-wave enter- 
ing the if* (892) selection. We model this component us- 
ing the LASS model [l?], EH which accounts for the inter- 
ference between the if,* (1430) resonance and the nonres- 
onant component. The shape of the Kir invariant mass 
is kept fixed to the results found in [17| : we fit for the 
LASS yield in the range 1035 < m(Kir) < 1550 MeV 
and extrapolate the expected yield to the signal region 
755 < m(Kir) < 1035 MeV. We find yields that are con- 
sistent with zero, ranging from -56 to 65 events. We fix 
this yield to zero if it is negative and take the estimated 
value otherwise. 

PDF shapes for signal, Kir, and BB backgrounds are 
determined from fits to MC samples, while for the qq 
background we use data samples from which the signal 
region, 5.27 < m ES < 5.29 GeV and \ AE\ < 0.075 GeV, 
is excluded. The calibration of m^s and AE is validated 
using high-statistics data control samples of B decays 
to charmed mesons with similar topologies (e.g. B — » 
D(Kirir)ir, B -> D(Kirir)p). 

We use linear combinations of polynomial, exponen- 
tial, and Gaussian functions to parameterize most of the 
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PDFs. For qq background, we adopt a parameterization 
motivated by phase-space arguments [19j |. 

We allow the most important parameters of the qq 
background to vary in the fit, along with the signal yield. 
Given that the signal yields we extract are small, we can- 
not vary the longitudinal polarization fraction //,. Since 
no strong theoretical predictions exist about its value, we 
impose /l — 0.5 and vary it within the physical range to 
evaluate the systematic uncertainty. We do not include 
the SXF component in fits with signal yields that are 
consistent with zero to avoid instabilities in the SXF fit- 
ted yield. In the case of the b\K*° mode, where the 
(statistical only) signal significance exceeds three stan- 
dard deviations, we retain the SXF component, fixing its 
yield to correspond to the rate given by the simulation 
for its size compared with the signal yield. In this case, 
introducing the SXF component causes the signal yield 
to vary by a small fraction of the statistical error. 

To evaluate the potential bias Y$ that arises from ne- 
glecting the correlations among the variables entering the 
fit, we perform fits to ensembles of simulated experi- 
ments. Each such experiment has the same number of 
signal and background events as the data; qq events are 
generated from the PDFs, while for the other categories 
events are taken from fully simulated MC samples. 

We compute the branching fraction B for each mode 
by subtracting Yq from the fitted signal yield Y and di- 
viding by the efficiency e and the number of B mesons in 
our data sample. We assume the branching fractions of 
the T(AS) to B+B~ and B°B° to be each 50%, consistent 



with measurements We evaluate e from signal MC 
samples, taking into account the difference in reconstruc- 
tion efficiency for longitudinally and transversely polar- 
ized events. For the K* + modes, we combine the branch- 
ing fraction results from the two sub-modes by adding 
their — 2 ln£ curves. The significance S is computed from 
the difference between the value of —2 In C at zero signal 
and its minimum value. The results are summarized in 
Table HTl while in Fig. [T]we show the projection plots onto 
the 771es variable for the ten final states we investigated. 
We do not observe a statistically significant signal for 
any of the eight decay modes. We quote upper limits on 
their branching fractions at the 90% C.L., taken as the 
branching fractions below which lie 90% of the totals of 
the likelihood integrals, constraining the branching frac- 
tions to be positive. The systematic uncertainties are 
taken into account by convolving the likelihood function 
with a Gaussian of width corresponding to the total sys- 
tematic uncertainties. 

We study the systematic uncertainties due to imper- 
fect modeling of the signal PDFs by varying the rele- 
vant parameters by their uncertainties, derived from the 
consistency of fits to data and control samples (the sys- 
tematic uncertainty on the signal yield varies from 0.6 
to 4.1 events, depending on the final state). The uncer- 
tainty due to the bias correction is taken as the sum in 



quadrature of half the correction itself and its statistical 
uncertainty (0.4-7.5 events). We vary the yield of the BB 
backgrounds by ±50% (the resulting uncertainty is 0.1- 
8.5 events) and the yield of the S- wave Kir component 
by the larger of ±100% of the extrapolated yield and its 
statistical uncertainty (0.2-14.3 events). The asymmet- 
ric uncertainty associated with fi, is estimated by taking 
the difference in the measured B between the nominal fit 
(/i = 0.5) and the maximum and minimum values found 
in the scan along the range [0, 1]. We divide these values 
by \/3, motivated by our assumption of a flat prior for 
/l in its physical range; this is one of the largest sources 




■.25 5.26 5.27 5.28 5.29 
m ES (GeV) 



25 5.26 5.27 5.28 5.29 
m ES (GeV) 



FIG. 1: Projections onto mES for the modes (a) b 1 p + , (b) 
b° lP +, (c) b+p°, (d) b° lP °, (e) b;K*+ +n0 , (f) blK*+ +v0 , (g) 

&r*£o*+' (h) b ° K K + ^+> 00 b t K *°> (j) b ° K *°- Points with 

error bars represent the data, the solid (dashed) line repre- 
sents the total (sum of the backgrounds) fitting function. The 
background is suppressed by a cut on ln£, optimized sepa- 
rately for each final state. 
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TABLE II: Signal yield Y and its statistical uncertainty, bias Yo, detection efficiency e, significance S (including systematic 
uncertainties) and central value of the branching fraction B with associated upper limit (U.L.) at 90% C.L. The efficiency e 
takes into account the product of the branching fractions of the intermediate resonances. 



Mode 


Y 


Yo 




S 


B 


B U.L. 




(evts) 


(evts) 


(%) 


{<?) 


(io- 6 ) 


(io- e ) 


b~p+ 


-33 ± 10 


4± 2 


3.0 




-1.8 ±0.5 ±1.0 


1.4 


6?P + 


-18 ±5 


-4 ±2 


1.1 




-3.0 ±0.9 ±1.8 


3.3 


6+p° 


37 ±25 


8±4 


3.6 


0.4 


1.5 ± 1.5 ± 2.2 


5.2 


b° lP ° 


-8 ±19 


5±3 


2.4 




-1.1 ±1.71^ 


3.4 



&-7T+ 1.7 2.4±i;| ± 1.0 5.0 

blK*+ +v0 3±8 -5 ±3 0.8 0.9 1.8 ±1.9 ± 1.4 

b^K*+ l + 17 ±9 4± 2 0.9 1.5 3.2 ±2.1±i'° 

b\K* + ^ 0.1 0.4t L °t 3 2 ° 6.7 

b°K*+ +v<3 -8 ±7 -3 ±2 0.5 - -2.2 ± 3.0±|;§ 

b\K*+ l + 3±4 0±0 0.4 0.4 1.6 ±2.5 ±3.3 

bfK*° S 55 ±21 15 ±8 2.8 1.5 2.9 ± 1.5 ± 1.5 5.9 

blK* 30 ±15 -6 ±3 1.7 2.0 4.8 ± 1. 9±H 8.0 



of systematic uncertainty, ranging from 0.1to3.6xl0~ 6 . 
Another large source of uncertainty is imperfect knowl- 
edge of the SXF fraction; based on studies of control 
samples performed in similar analyses, we assign a 5% 
multiplicative systematic uncertainty on the SXF frac- 
tion (relative to correctly reconstructed signal) for each 
7T° in the final state. Other uncertainties arise from the 
reconstruction of charged particles (0.4% per track), K® 
(1.5%), and 7r° mesons (3% for 7r°); the uncertainty in 
the number of B mesons is 1.1%. 

In summary, we present a search for decays of B 
mesons to b\p and b\K* final states. We find no sig- 
nificant signals and determine upper limits at 90% C.L. 
between 1.4 and 8.0 x 10~ 6 , including systematic uncer- 
tainties. Though these results are in agreement with the 
small predictions from naive factorization calculations 
, they are much smaller than the predictions from the 
more complete QCD factorization calculations 0- The 
fact that the branching fractions for these AV modes are 
smaller than our previously measured AP modes is 
surprising given that the opposite is expected based on 
the ratio of the vector and pseudoscalar decay constants. 
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